Enhanced rise of rogue waves in slant wave groups 
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Numerical simulations of fully nonlinear equations of motion for long-crested waves at deep water 
demonstrate that in elongate wave groups the formation of extreme waves occurs most intensively if 
in an initial state the wave fronts are oriented obliquely to the direction of the group. An "optimal" 
angle, resulting in the highest rogue waves, depends on initial wave amplitude and group width, and 
it is about 18-28 degrees in a practically important range of parameters. 
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The phenomenon of rogue waves (also known as freak, 
killer, giant, or extreme waves) at the ocean surface has 
attracted much attention in recent years (an extensive 
discussion can be found in [H-ll] , and particular aspects of 
the rogue wave formation are considered, for instance, in 
[i-[l3). A probable scenario suggests that linear mecha- 
nisms, such as interaction of surface waves with a nonuni- 
form current, cause preliminary amplification of wave 
amplitude, making the most wide and tall wave groups 
unstable with respect to the so called modulational insta- 
bility (20l - [2^ . A rogue wave is thus the final stage of the 
development of that instability, as it has been confirmed 
by direct numerical simulations of exact equations of mo- 
tion for potential two-dimensional (2D) flows of a perfect 
fluid with a free surface [23.-.26j . However, many three- 
dimensional (3D) aspects of the problem still remain far 
from being clear. Partly the difficulty is explained by 
the absence of compact and explicit exact equations of 
motion for a 2D free surface in the 3D space, which fact 
results in rather slow and cumbersome implementations 
of the existing numerical methods based on the Euler 
equations. Therefore some approximate analytical and 
numerical models were suggested to study 3D dynam- 
ics of oceanic waves. In particular, for weakly nonlinear 
regimes the nonlinear Schroedinger equation (NLSE) 
and its generalizations [l^ [H, i27 29] are widely used, 
which are simplifications of the Zakharov equation tak- 
ing into account renormalized 2—^2 wave processes, 
which equation is the most general among weakly non- 
linear models (see [2l|, H^l , and references therein) . But 
weakly nonlinear equations are definitely not appropriate 
to describe rogue waves at their final stage of evolution. 
That is why another, fully nonlinear approximate model 
has been developed by the present author, based on a 
different small geometrical parameter, which is a small- 
ness of deviation from a planar flow [3l|. In other words, 
a narrow angular distribution of the wave spectrum in 
the horizontal Fourier plane is assumed. The model is 
good for arbitrary steep long-crested water waves propa- 
gating closely to x direction in the horizontal (x, q)-plane 
{y is the vertical coordinate) . The corresponding numer- 



ical method is reasonably fast [32], and it was applied 
to study breathing rogue waves in a random wave field 
[ss'] , nonlinear stage of the modulational instability with 
specific zigzag coherent structures producing freak waves 
through mutual interactions [34[, and the two kinds of 
rogue waves in weakly-crossing sea states [s^ . Additional 
numerical examples can be found in the recent paper [36| . 

In the present work the author continues investigation 
of 3D effects in the dynamics of extreme water waves. 
Let us consider sea states in situations when a typical 
wave has a length A and an amplitude A, so a typical 
wave steepness is s = 2'kA/\. As a general rule, a width 
of the wave group NX should be sufficiently large for the 
rogue wave formation to take place. 
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with some constant /c ~ 1 (precise value of is not im- 
portant, since we speak now about typical values; besides 
that, there exist several slightly different definitions for 
rogue waves; see discussion in [4!|). The product sN is the 
so called Benjamin-Feir Index (BFI) With small BFI, 
dispersive effects in the wave dynamics are dominant, 
while for large BFI the nonlinearity becomes essential. 
The above condition comes from simplified consideration 
of the phenomenon within one-dimensional (ID) NLSE 
describing a complex wave envelope B{x,t) in the case 
of a planar flow (see [2ll|). 
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where k = 2t:/X is the wave number, w = vy^ is the wave 
frequency, and g is the gravitational acceleration. Con- 
dition (m means that the wave group contains in some 
sense a soliton of the ID NLSE. By the way, the pa- 
rameter TV is the number of individual waves in the ID 
group. One of the most intriguing questions in the the- 
ory of 3D rogue waves is about their appearance in non- 
coherent random sea states, when typical wave groups 
are not very tall and/or wide, so sN < Ic- However, 
in the cases when envelopes B{x,q,t) of wave groups 
have a length much longer than the width (elongate wave 
groups, as in weakly-crossing sea states [ssj; see also 
[3 7 1] ) , then an additional important parameter comes into 
play, namely an angle 9, at which wave fronts are ori- 
ented relatively to the "long" direction. The purpose of 
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the present work is to demonstrate via numerical sim- 
ulations that in this regime the most high freak waves 
arise for some optimal angle 9, which depends on the 
parameters N and s. This essentially three-dimensional 
effect is reported for the first time, and it is very dis- 
tinct in practically important ranges 0.10 s < 0.14 
and 5 ^ N 10. In this parametric region we find 
9 — 18 • • • 28°, and arising rogue waves have the height 
^max ~ 0.06A at which the process of wave breaking be- 
gins to occur, while the ratio FmaxM = (2.5 • • -3). At 
the same time, ymin ~ — 0.04A, due to the crest-trough 
asymmetry of gravity water waves. 

To deal with a minimal set of parameters in our study, 
we consider idealized wave groups which are infinitely 
long in one direction. So we have stripes along X2 
axis in a turned {xi,X2) coordinate system which is ori- 
ented at some angle 7 to the (x, q)-system, with 7 be- 
ing slightly less than 9. The components of the corre- 
sponding wave vector in (xi,X2) coordinate system are 
(fc cos 0, — fc sin 0). Thus, initial wave fronts are oriented 
not exactly parallel to g-axis, but at a small angle (about 
several degrees) clockwise, while the stripe itself is ori- 
ented at the angle 7 anticlockwise. This is made because 
crests of arising extreme waves are always oriented more 
gently to the stripe direction comparatively to the crests 
in the initial state (see Fig.l, and also [SJ, H^), so the 
choice 7 < results in a more close orientation of rogue 
wave crests to g-axis, as it is required for applicability of 
the employed approximate quasi-2D model [3ll |32||. 

In the initial state, the complex envelope of the first 
wave harmonics is put purely real and given by a simple 
expression, 

4-^(2^)- 

Thus, we can identify the parameter N as follows: N w 
4w. We also add a low- level random-phase perturbation 
into initial wave spectrum, similarly to [ssl Is^ j. to be 
sure that our results are robust with respect to a noise 
in initial conditions. 

For convenience of graphical presentation we choose 
A — 100 m, so the corresponding wave period is T = 
y^2TrX/g 8 s. The computational domain has the 
rectangular shape Lx x Lq, with the periodic bound- 
ary conditions in both directions, and = 2 km. For 
the parameter Lq, several different values were taken 
{Lq = 4, 5, 6, 7, 8 km) in order to ensure the quasi-2D 
regime for different angles 0, at least in an initial stage 
of evolution; by the way, 7 — a,Tcta,ii{Lx/ Lq). Exception 
is for small 9, when Lq = Lx, and 7 = 0. For example, 
in Fig.l shown are two sub-regions of the whole domain 
2 km X 7 km. 

The simulations were performed on modern personal 
computers using numerical method described in [32| . The 
final resolution was about 12000 x 3000 points in the cases 
when extreme waves evolved closely to breaking (the be- 
ginning of the breaking is characterized by a rapid in- 
crease of the maximal wave steepness after reaching a 
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FIG. 1: Numerical example of formation of rogue waves in a 
slant wave group. 



critical value s* ~ 0.5 rad, which is about the steep- 
ness of the limiting Stokes wave). Some of the obtained 
numerical results are presented in the figures. In partic- 
ular. Fig. la shows a map of the free surface &t t — Q, 
for s = 0.14, TV w 6, and w 18.4°, while Fig.lb is a 
map for a later time moment, about several tens of wave 
periods, when rogue waves and deep troughs form a spe- 
cific slant structure resembling wake waves after a ship. 
It is worth noting that fragments of similar wave stripes 
develop spontaneously in nonlinear stage of the modula- 
tional instability, where they form zigzag patterns, with 
rogue waves arising mainly at zigzag turns [13, Ha] . In 
Fig. 2 presented are some wave profiles from Fig.l, which 
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FIG. 2: Some wave profiles corresponding to Fig.l. Tlie rogue 
waves concentrate tlie energy, and therefore amplitude of the 
remaining waves in the group is decreased. 



emphasize amazing features of rogue waves, such as their 
strong localization and a high relative amplification. In 
general, a highly nonlinear nature of the rogue wave phe- 
nomenon is confirmed in our numerical experiments. 

Very interesting and important is Fig. 3, where depen- 
dences of the maximum surface elevation versus time for 
different 9 are compared. It is clearly seen that for s and 
w satisfying the relation Asw « 1 (nearly critical BFI), 
the most "perfect" rogue waves with Fmax ~ 0.06A arise 
if 9 is sufficiently large, about 18 • • • 28°, while for small 
6 just a moderate wave growth takes place, with subse- 
quent decrease. It should be noted here that in the ab- 
sence of wave breaking, after the decrease, the second rise 
is observed at later times (not shown), but in this work 
we do not discuss that quasi-recurrent behavior caused 
by complex interaction of quasi-solitons constituting the 
wave group. 

An optimal angle 9{s,w) increases as s or u; decreases 
(compare Fig. 3a to 3b, and 3c to 3d). However, there ex- 
ists a critical value 0* = arctan(l/\/2) ~ 35.3°, at which 
the second-order dispersive coefficient of the correspond- 
ing ID NLSE, approximately describing the dynamics of 
the wave envelope B{xi,t) in a moving frame of refer- 
ence, changes the sign (see [HI, HH] for details): 
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Our simulations show that actually for 9 approaching 
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FIG. 3: Maximum elevation of the free surface versus time in 
slant wave groups with fixed width and steepness, for different 
angles Q. Waves shown in Fig. lb and Fig. 2b correspond to 
dependence '3' in Fig. 3c, slightly before the end point. Wave 
breaking takes place at the end of experiments a5, b3, b4, b5, 
c3, c4, c5, c6, d3, d4, and d5. 
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the most tall waves in the stripe become significantly 
shorter than A, and finally they break with l^ax well 
below the value 0.06A. It should be also noted that for 9 
close to 6^, the quasi-2D regime is violated after a short 
time (crest orientation is then strongly varied across the 
stripe; not shown), and therefore that parametric region 
cannot be accurately investigated with the help of the 
quasi-2D model. 

Qualitatively, such enhanced growth of extreme waves 
in slant wave groups can be explained by noting that with 
9 one should modify the condition ([1]) for rogue wave 
formation as follows, 

sN/\/cos^9-2sm^9 > Ic, (5) 

since the spatial scale perpendicular to the stripe is for- 
mally renormalized in the ID NLSE by the above square 
root, as it is clear from Eq.Q. However, Eq.(|31) is not 
applicable with small values of the root, because the nec- 
essary condition of spectrally narrow wave field is then 
violated very soon in the course of evolution. Adding 
higher-order linear dispersive terms into Eq.Q, one can- 
not improve situation, since the rapid widening of wave 
spectrum is a real physical effect for 9 close to 6',, and 



thus the nonlinear term in NLSE should as well be modi- 
fied to a non-local form determined by the 4- wave matrix 
element of Zakharov equation [111, [s^. We do not write 
here the corresponding ID reduction of the Zakharov 
equation for slant wave stripes, since it is too difficult 
for analytical treatment, and besides that, it is not ac- 
curate for large wave amplitudes. So, in the absence of 
reliable analytical estimates for dependence of the max- 
imum surface elevation on s, N, and 9, the above pre- 
sented numerical results are quite valuable. 

To summarize, it has been shown in this work for the 
first time that an oblique orientation of wave fronts in ob- 
long water-wave groups is able to enhance the process of 
formation of rogue waves. This 3D effect is most promi- 
nent when the Benjamin- Feir Index of the wave field is 
slightly below its critical value. 
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